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O ■ ABSTRACT 

^ ■ 

Aims. In a previous work, Ulysses data was analyzed to build a complete axisymmetric MHD solution for the solar wind 
at minimum including rotation and the initial flaring of the solar wind in the low corona. This model has some problems in 
reproducing the values of magnetic held at 1 AU despite the correct values of the velocity. Here, we intend to extend the previous 
analysis to another type of solutions and to improve our modelling of the wind from the solar surface to 1 AU. 
O I Methods. We compare the previous results to those obtained with a fully helicoidal model and construct a hybrid model combining 
both previous solutions, keeping the flexibility of the parent models in the appropriate domain. From the solar surface to the 
Alfven point, a three component solution for velocity and magnetic field is used, reproducing the complex wind geometry and 
the well-known flaring of the field lines observed in coronal holes. From the Alfven radius to 1 AU and further, the hybrid model 
, keeps the latitudinal dependences as flexible as possible, in order to deal with the sharp variations near the equator and we use 
• • • the helicoidal solution, turning the poloidal streamlines into radial ones. 
_ ' Results. Despite the absence of the initial flaring, the helicoidal model and the first hybrid solution suffer from the same low values 
, of the magnetic field at 1 AU. However, by adjusting the parameters with a second hybrid solution, we are able to reproduce 
» I ■ both the velocity and magnetic profiles observed by Ulysses and a reasonable description of the low corona, provided that a 
certain amount of energy deposit exists along the flow. 

Conclusions. The present paper shows that analytical axisymmetric solutions can be constructed to reproduce the solar structure 
and dynamics from 1 solar radius up to 1 AU. 
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1. Introduction 



Since iPar keif s model (1958), many studies have been pre- 
sented to explain and predict the features and properties 
of the solar wind, mainly following two different, yet com- 
plementary, approaches, kinetic and fluid approximations. 

These techniques are able to reproduce certain aspects 
of the observed solar wind but both show some limita- 
tions, mainly due to the complexity of the several acceler- 
ation mechanisms, the uncertainties concerning the origin 
of the fast solar wind, the associated problem of coro- 
nal heating, etc. Different models have been presented 
improving results of the acceleration. Tw o fluid m odels 
and more recently three-fluid models, (e.g. lOfmanl . 2000; 
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Zouganelis et al.l . 2004) have been constructed to explore 
the kinetic aspects of the wind acceleration by supra- 
thermal electrons in the collisionless region far from the 
Sun. All models still have difflculties avoiding very high 
temperature for the electrons. Other s o urces of heating 
such a s turbulence (Landi fc PanteUini . 20031) or Alfven 



waves (jUsmanov fc Goldsteinl . l2003l ; iGrappin et al.l I. l2002h 
may also explain the acceleration by lowering the effective 
polytropic index of the flow. This point is not yet resolved 
and we shall not address this question here. Instead we will 
invoke the need for turbulence or Alfven wave damping in 
our solutions. 

Another approach consists of constructing MHD so- 
lutions to analyze the 3D structure of the wind, al- 
most independently of the heating source. Various mod- 
els have been constructed, either 2-D ones able to de- 
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scrib e the low corona o f the s un up to 10 sol ar 
'e.g. Pneuman &: Kopnl 1971; Steinolfson et al 



radii 

^ - . . - - , 1982; 

Cuperman et all 1990 and references therein), the slow 
solar wind inside the brlRhtness b oundary 



(e.g. Nernev fc Suessl 2005 and refer- 



nal streamers, 
ences therein), 2-D ones for larger distances and mod- 
els for all range of distances. Those models were pro- 
posed because the flaring of the streamlines in polytropic 
winds favors the acceleration. Recent observations by 



Wang fc SheelevI (|2003l ) showed, however, that this may 
not be the case for the real solar wind. This favors a de- 
scription of the 3D structure of the solar wind using self- 
similar MHD analytical solutions for non poly tropic winds 



(|Tsinganos fc Sautvl . 119921 : iLima et al.l . l200lh . In the first 



of these two models it has been shown that the flaring of 
the lines may instead limit the acceleration of the wind. 

An increasing amount of observational data is now 
available. Ulysses measured for the first time the mag- 
netic field, the dynamics and the temper ature of the wind 



around 1 AU out of the ecl iptic plane (IMcComas et al 



2000l). Data from ACE (jStone et all . Il998l). WIND 



(|Acuna et all [l995|) SoHQ UVCS (e .g. IWoo fc Habbll . 
20051 LASCO fe.g. lLewis fc Simnettl . 2002) and SoHO 
CDS (jOallagher et all . Il999h are providing new insights 
into the origin of the solar wind within coronal holes . 
Doppler Scintillation measurements ( Woo fc Gazid . 1994D 
also brought new constraints to solar wind modelling. 

Semi-empirical models that use data to set bound- 
ary conditions for a numerical approach to the problem 



have also been propose d (e.g 



Sittler fc Guhathakurtal (|l999l ): ICroth et al 



Steinolfson et al.l ([19821); 

0999)). 



Nevertheless some doubts o n the boundarie s of so me sim- 
ulations are still present ( Vlahakis et al. . 2000l ). More 
recently, some new developments have suggested that 
numerical simulations ca n benefit greatly from an an- 
alytical treatment (e .g. iKep pens fc Goedbloed . 2000; 
Usmanov fc Goldsteinl . 2003; Jlayashi, 2005). Numerical 
simulations are still quite time-consuming although this 
is rapidly improving. However, there are other limitations 
such as maintaining divergence-free magnetic fields, limit- 
ing the effects of numerical magnetic diffusivity or solving 
the 3D structure of the wind including rotation even at 
large distances. Note that the main problem with present 
simulation s is the existenc e of a numerical magnetic diffu- 
sivity (e.g. Grappin et al.l . 2002). This is why we propose 
to construct semi-analytical models which are less sophis- 
ticated than numerical simulations but simpler to handle 
and more versatile. They also provide a complementary 
approach. 

In the p resent work, that follows closely the work of 
Sautv et "all (2005 hereafter SL1T05) we focus on the dy- 
namics of the protons in the solar wind. We apply known 
MHD analytical models to Ulysses data at solar mini- 
mum and test their advantages and limitations. Wc will 
generate an exact wind solution based on the model of 
Lima et al. (2001 hereafter LPTOl) and on the data fit 
made in SL IT05. In SLIT05 two models (LPTOl and 
Sautv et all 1999 hereafter STT99) were used to fit the 



data. The final solution was based on solving the differen- 
tial equations of the latter of these two models. Finding 
a solution that complies with the constraints given by the 
data fit and the ones from its topological features is not 
easy. 

Regarding the limitations of LPTOl wind solution 
(purely radial, yet very adaptable to the latitudinal de- 
pendences) and the ones of the similar study presented 
in SL1T05 (see also ST T99; ISautv et all 2002, hereafter 
STT02 and [Sautters, 2004, hereafter STT04), we take 



into account the advantages of both models by creating 
hybrid solutions. These use the 2.5D features of the STT04 
model to describe the solar wind dynamics from the solar 
surface towards the Alfven sphere. From the Alfven point 
towards 1 AU and beyond, these hybrid solutions will use 
the advantages of the LPTOl model in fitting steep varia- 
tions of velocity, density and magnetic field with latitude 
and expressing the radial behaviour of the solar wind in 
this region. However, we still solve the complete set of 
MHD equations in the radial domain and not simply the 
Bernoulli equation along the streamlines. Thus, the solu- 
tion remains consistent everywhere. We discuss the prop- 
erties of the solutions thus obtained and physical grounds 
for their limitations. 

We maintain the criteria used to find a good solution 
from the ULYSSES data fit and the measured values of the 
physical quantities at 1 AU. It will be shown later that, for 
some sets of parameters, both the LPTOl model and the 
first hybrid model show the same problems mentioned in 
SLIT05, namely in reproducing the values of the magnetic 
field at 1 AU from Ulysses. These will be solved by a 
judicious choice of parameters in the second hybrid model 
that generates a solution consistent with Ulysses data. 



2. Self-similar MHD outflow models from central 
rotating objects 

The following two axisymmetric wind models are obtained 
by self-consistently solving the full system of ideal MHD 
equations. In the present work we use spherical coordi- 
nates [r, 6, (j)]. All quantities have been normalized to their 
values at the Alfven radius along the polar axis, similarly 
to SLIT05. They will be identified by the subscript *, i.e. 

, and for velocity, density and magnetic fields at 
the Alfven polar point, respectively. All equations will be 
presented in a normalized form where the distance to the 
solar surface is related to the real distance by i? = r/r*. 

2.1. Model A with flaring streamlines 

In model A all three components of the velocity and mag- 
netic fields are accounted for (STT99, STT02, STT04). 
Nevertheless, an expansion up to first order in latitude 
of the forces is performed by using harmonics with polar 
values as references. Such a procedure makes the whole 
system analytical tractable and also describes the helio- 
latitudinal variations of the wind quantities. The fields 
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describing the outflow dynamics are 

Ve{R,9) = -K 
V^{R,9) = 
Bj. (^R, 9 
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(3) 

(4) 
(5) 
(6) 
(7) 
(8) 



where Vr, Ve, are the three components of the velocity 
field, Br, Bg, B^, the three components of the magnetic 
field, p, the density, P, the pressure and C is a constant. 
There are three functions of i?, namely M, / and 11. 

2.2. Model B with helicoidal/ radial streamlines 

Model B assumes a simpler geometry with radial stream 
and field lines in the poloidal plane {i.e. zero 9 components 
of the velocity and magnetic fields). It is more versatile at 
reproducing steep latitudinal variations (LPTOl). In this 
case, the fields describing the outflow dynamics are 



2e , 



Vr{R,9 

V^{R,9 

Bj. (i?, 9 
B^{R, 9 
p{R,9 
P{R,9 



1 + fJ. sin' 
1~AP/R^\ Rsin^e 
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I + 68111^' 9) 



Rsm'9 
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^P*K'(no + ni sin^^i 



(9) 
(10) 

(11) 

(12) 

(13) 
(14) 



where the same notation is used and Ho and Hi are func- 
tions of R. 



2.3. Geometry of the solutions 

The relevant wind type solutions cross various critical 
points related to the non-linearity of the system of equa- 
tions and its mixed elli ptic/h yperbolic nature (see for in- 
stance Tsinganos et al.. ll996t STT04). Each model is de- 
scribed by three functions of R, M(R), Il{R) and f(R) for 
model A, M(R), Uo{R) and Ili{R) for model B. The func- 
tion f(R) characterizes the geometry of the fieldlines and 



Table 1. Parameters obtained by data fitting at 1 AU 
along the polar axis using the Ulysses hourly average 
data, except for -Bt, lAU which is calculated along the equa- 
torial plane. 



Parameters 


Model B 


k/iau 


0.35 


<5/lAU 


1.95 


/lAU 


1.00 


£ 


5.64 




- 0.18 


ViAu (km/s) 


775 


niAu (cm^^) 


2.48 


BiAu (m G) 


30.4 


Bt,iau (m G) 


29.5 



expresses the expansion factor. For a fully radial poloidal 
fieldline (i.e. an helicoidal pattern of the lines in 3D) we 
have f = 1 which is the case of model B. The func- 
tion M{R) describes the poloidal Alfven Mach number 
which is unity at the Alfven radius. At this point the ki- 
netic energy overtakes the magnetic one. A limitation of 
both models comes from their self-similar nature. Thus the 
Alfven Mach number is independent of latitude and there- 
fore the Alfven iso-surfaces are spherical. The functions 
n, Ho and Hi are determined by numerical and analytical 
techniques that are explained in STT02 and LPTOl. 

3. A complete solution with helicoidal/radial 
streamlines 

In SLIT05 we fitted Ulysses data using the latitude de- 
pendence of models A and B. We have shown that both 
models yield similar parameters. The system of ODEs was 
integrated exclusively using model A. Conversely, in this 
section we use the ODEs of model B to derive a full solu- 
tion and compare the results with the ones from SLIT05. 
The model flexibility provides a better fit of the latitu- 
dinal functions, which may be crucial in dealing with the 
poloidal data at 1 AU. Yet, this solution cannot reproduce 
the flaring of the streamlines as they remain radial in the 
poloidal plane. 

3.1. Method for a solution 

The free parameters of the model (e, 8 and /i), the polar 
values of the number density, radial velocity and magnetic 
field at 1 AU, niAu, ^lAU and BiaUj respectively, and the 
equatorial toroidal magnetic field at that same distance, 
-Bt,iaUj have already been constrained by the Ulysses 
data fitting procedures used in SLIT05. The end results 
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are summarized in Table [TJ From these constrained pa- 
rameters the latitude functions are well defined. The radial 
dependences of the physical quantities are determined by 
integrating the ODEs of the model. From the knowledge 
of the values at 1 AU of the polar velocity, Viau , density, 
piAU, and magnetic field, i?iAU, it is possible to infer the 
value of the Alfven number at that same distance, 



47rpiAulfAU 



r2 

-°1AU 



Simultaneously using Eq. 
tant reference value. 



(15) 



T5)) we obtain another impor- 



p* 



lAU * lAU 



-°1AU 



(16) 



In the original paper (LPTOl) the relations that rule 
the model arc normalized to the solar surface. The critical 
solution for the solar wind is calculated based on three 
simple criteria. The first one is Vr{r/rQ = 1)/Vq = 1 and 
the other two are the continuity in the acceleration at the 
Alfven singularity and fast magnetosonic separatrix. In the 
present work we keep the two last criteria but Eqs. ([9]) to 
are normalized to the Alfven radius as in SLIT05. 
Thus, regarding the first criterion, the definition of the 
solar surface, r^jr^ = Rq, poses a problem. In order to 
calculate it we used the value at which the radial velocity 
goes to zero or reaches its minimum value. As we expect 
the radial velocity to have a very steep variation at those 
distances, the corresponding error at the evaluation of the 
solar surface will be very small. The final solution should 
also reproduce the measured values obtained by Ulysses 
at 1 AU (mainly radial velocity, radial and toroidal mag- 
netic field and density). The total acceleration between 
the Alfven radius and 1 AU can be parameterized by, 

^lAU 



V 



(17) 



Guessing this parameter, we obtain an initial value of U* 
and -B* = y/iiFp^V^, . From Eq. ([TT|) it is possible to deter- 
mine the value of 1 AU in Alfven radius units, 



RlAV 



lAU 
^(AU) 



215r0 

''*(r0) 



B 



lAU 



(18) 



The last equation gives the location of the solar surface 
since -Riau — t'iau/''* — 215 i?©. Already having the 
anisotropy parameters, e, d and /i, we still need A and 
V. Combining Eqs. pT|) and (H^, assuming that we are 
at large distances (i?iAU >> 1) and since the lines are 
radial, which is very reasonable at 1 AU, we get, from Eq. 
(fT2|) applied at 1 AU, on the equatorial plane. 



A 



B' 



T.IAU 



lAU 



i?lAU 

By definition, we also have 



l2GM 



(19) 



(20) 



Table 2. Input and output data for the critical wind so- 
lution calculated using model B, hybrid 1 and hybrid 2. 





input param. 


model B 


hybrid 1 


hybrid 2 




r 



l.yo 


o on 






e 


5.64 


5.64 


5.64 




P 


- 0.682 


-0.406 


-0.029 




K 




0.20 


0.0123 


At 




2.15 


2.80 


1.90 


1 AU 


X 


0.1662 


0.2468 


0.1383 






1.462 


0.8872 


0.3767 




ViAu (km/s) 


775 


775 


775 




niAu (cm~^) 


2.48 


2.48 


2.48 




BiAu (^G) 


2.432 


9.81 


30.4 




-Bt.iau ifJ-G) 


29.5 


29.5 


29.5 


output param. 


At 


V = ViAv/V, 


2.151 


2.193 


1.908 


1 AU 


RlAV 


156.9 


137.0 


13.37 




ViAu (km/s) 


775.4 


607.6 


779.4 




niAu (cm"^) 


2.413 


0.207 


2.414 




BiAu (pG) 


2.369 


0.641 


29.8 




-Bt.iau (mG) 


29.45 


9.64 


29.1 




TiAu (10" K) 


3.037 


4.246 


41.88 


At the 


V, (km/s) 


360.5 


277.1 


360.5 


Alfven 


n. (10^ cm~^) 


131.3 


8.095 


0.8423 


radius 


B, (10^ /iG) 


5.989 


1.143 


0.5435 




St* (10^ fiG) 


9.953 


1.360 


0.278 




T. (10® K) 


9.892 


3.902 


3.597 




Fig. 1. Poloidal fieldlines and density contours of the so- 
lution of Table [2] for model B. Distances are given in solar 
radii. The solid circle line indicates the Alfven singularity 
and the fast magnetosonic separatrix which are almost 
coincident. 
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1 10 100 1.3681 1.3691 1.3701 1.3711 

distance (solar radius) distance {soiar radius} 



Fig. 2. Panel (a) polar radial velocity as a function of distance. Panel (b) a zoom view of the singularities: the 
Alfvenand the fast critical points, represented by * and x respectively. 



Our procedure for finding a critical solution, using 
model B, that fulfills all criteria is thus: 

- By fitting ULYSSES data at 1 AU obtain the anisotropy 
parameters, e, 5 and /.t, and the values BiaUj ^t.iaUj 
PiAu and ViAu; 

- calculate Rq, A and v with /iau 1 and using Eqs. 

(HHD to (HOI); 

- with an initial guess of -q, determine the values of 14 
and from Eqs. ^ and (fTe]) : 

- at this stage it is possible to build a critical solution us- 
ing the criteria of acceleration continuity at the critical 
points; 

- this solution will give new values for the solar surface 
radius and for the different physical quantities at 1 AU; 

- iterate until the computed values of the solar surface 
radius and velocity at 1 AU are close to the fitted ones. 

Two convergence criteria are inherent to this proce- 
dure, the distance of the Alfven surface above the Sun (or, 
equivalently, the value of the magnetic field strength at 1 
AU - see Eq. [T8|) and the velocity at 1 AU. Satisfying all 
the criteria only by changing 77 is not possible. Therefore, 
this can only be achieved by changing, in addition, at 
least one of the parameters, thus releasing one of the con- 
straints. Considering that there are five constraints given 
from the data fit, exploring all the parameter space is a 
formidable task. The set of parameters concerning the best 
possible solution is presented in Table [21 This will be dis- 
cussed in the following section. 

3.2. Results 

As can be seen in Table [2] we have constructed a solution 
where the convergence criteria are quite well fulfilled with 
input/output ratios very close to unity for all physical 
quantities. In Fig. [1] we show the fieldlines and the density 
contours in the poloidal plane. In Fig. [2l where the verti- 



cal dashed line represents the Alfven radius, the accelera- 
tion at the critical points is clearly continuous. Note the 
presence of two different critical points, very close to one 
another in Fig. [2] (b). However, searching for a converged 
solution led us to this single set of parameters by using a 
value of the radial magnetic field, -BiaUi of the order of 
one tenth of the value measured by Ulysses. A similar 
discrepancy was also found in SLIT05 using the STT99 
model instead of LPTOl model. Moreover, the parameter 
that evaluates the anisotropy of the radial magnetic field, 
fi, has also suffered a shift in its value (compare Table [T] 
to Table [2]) . As mentioned above we have tried to change 
the other input parameters and calculate their infiuence 
on the solution. The best option was to change those two 
parameters. Although these trials have been nearly ex- 
haustive, degenerated sets of the input parameters for the 
same critical solution may be possible. In Fig. [3] we show 
the temperature profiles at various latitudes. The higher 
effective temperature along the polar axis corresponds to 
the fast solar wind. At lower latitudes the lower temper- 
ature is related to a mixing between the fast and slow 
wind, which also corresponds to the lower velocities seen 
in Fig. [2] (a). The temperature distribution is similar to 
the one presented in SLIT05 although its maximum value 
is slightly better, around 10 x 10^. 

We have also calculated the effective polytropic index 
of this solution. After the temperature peak, it is almost 
constant, between 1.1 and 1.3 (see Fig. HJ. This v alue is 
quite close to the value inferred by Kopp and Holzer ( 19761 ) 
in their early polytropic model. Despite the difference be- 
tween model A and B in their geometry, neither can repro- 
duce from the observations the high value of the magnetic 
field inferred at 1 AU. The calculations were made such 
that we keep the temperature as low as possible and a 
reasonable value of the magnetic dipolar field at the base 
of the corona. We reproduce successfully the temperature 
and the velocity profile at 1 AU. It seems that the geom- 
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distance (solar radius) 



Fig. 3. Effective temperature as a function of distance for 
different values of latitude for the model B solution. The 
vertical dashed line represents the Alfvcn radius. 



10 100 
distance (solar radius) 



Alfven radius, and the other from the Alfven radius to- 
wards infinity. In the first domain we use model A, since 
the 3D structure is more able to describe the wind struc- 
ture near the solar surface. Although the latitudinal func- 
tions of model A are not very versatile it does provide a full 
3D description of the flaring. In the domain further out, 
where the fieldlines are almost radial in the poloidal plane, 
we can use model B and take advantage of its flexibility in 
fitting very steep variations of the physical quantities with 
latitude. The border between these two different domains 
was arbitrarily set at the Alfven surface. 

The major drawback of this construction is that we 
cannot guarantee continuity of all physical quantities ev- 
erywhere except along the polar axis. Generating the crit- 
ical solution with model A means that it must cross both 
slow magnetosonic separatrix critical point and Alfvenic 
singular point. In addition, the critical solution with model 
B has to cross the Alfven point and a fast separatrix crit- 
ical point. Thus, a new feature of this hybrid model com- 
pared to our previous solutions (SLIT05 and Sect. [3]) is 
the crossing of the three usual MHD critical points. Such 
a situation was present only in the over-pressured solu- 
tions presented in STT04. Since for the solar wind the fast 
point is close to the Alfven one, this was one more argu- 
ment to construct the hybrid solution starting precisely at 
this Alfvenic transition. Moreover, in order to match the 
two solutions we must search for continuity of the physi- 
cal quantities at the boundary as much as possible. This 
means that at the Alfven surface we ask for continuity of 
the density, pressure, velocity and magnetic field plus con- 
tinuity of the acceleration and fieldline geometry. Strictly 
speaking, this can only be done along the polar axis be- 
cause the latitudinal dependences of the physical quanti- 
ties are not identical in both models. Mathematically, we 
have 



model A 



odelB 



n 



jnodelA 



c = n; 



rnodelB 



(21) 
(22) 



Fig. 4. Polar polytropic index as a function of distance for 
the model B solution. The vertical dashed line represents 
the Alfven radius. 

etry does not control the decrease of the magnetic field 
but rather the temperature profile. It is even more sur- 
prising that in this LPTOl solution the density at 1 AU 
remains at a reasonable level, thus both the velocity and 
the mass flux at 1 AU correspond to the observed val- 
ues. It is thus more consistent to build an hybrid solution 
combining both models A and B. 

4. Hybrid solutions 

4.1. Method for a solution 

For the construction of the hybrid model we consider 
two different domains, one from the solar surface to the 



modelA 



modelB 



B 



modelA 



B 



modelB 



modelA mode 



IB 



dR\, 



dR\r, 



dR\ 



model A 



dR\ 



odelB 



(23) 
(24) 

(25) 
(26) 



The technique used to obtain a full hybrid solution is 
as follows. First, the radial velocity at the Alfven point, 

is determined by the same procedure as in the pre- 
vious section. Then, the model A critical solution has to 
cross the slow separatrix critical point and the poloidal 
fieldlines have to be radial at the Alfven point, Eq. (P^ . 
This yields the value of the velocity slope (i.e. the acceler- 
ation) at that transition point. Similarly, the value of the 
acceleration at the Alfven point is determined by crossing 
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Equeinriel ails 

Fig. 5. Plot of the fieldlines concerning the hybrid solution 
and density contours presented in Table[21 column 3, [hy- 
brid 1). Distances are given in solar radii. The black cir- 
cles corresponding to the slow magnetosonic critical sep- 
aratrix, the Alfven singularity and the fast magnetosonic 
critical separatrix, which almost coincide with each other. 



the fast separatrix critical point for the critical solution 
of model B. In order for this slope of the velocity to be 
equal on both sides of the transition point, Eq. (j25p . the 
value of /I had to be changed from the one derived from 
Ulysses data. However this parameter is not very well de- 
termined and affects only model B. In model A its value 
is fixed to —1 by construction and cannot be fitted. For 
model A, knowing the slope of the velocity, Eq. and 
the geometry, Eq. (|26p. at the Alfven point, the value of 
is fixed by crossing the Alfven point. Simultaneously, 
the value of XIq is fixed by the condition that the pressure 
is zero at infinity. This determines the constant C by Eq. 
(|22|l . Finally, the Alfven distance is fixed by the magnetic 
field strength at 1 AU for both models, Eq. (HH]). Thus, the 
value of K is determined such that the solution of model 
A matches the solar surface, Rq. 

4.2. Results for a hybrid solution - hybrid 1 

Table [2] shows the input and output values for the most 
important physical quantities regarding the first hybrid 
solution obtained using the technique presented in the pre- 
vious section. In this case convergence between the input 
and output parameters is less satisfying than in Sect. [3] 
This solution is hereafter referred to solution hybrid 1. 
Figure [5] shows the geometry of the fieldlines for this so- 
lution. It clearly shows that beyond the Alfven point the 
fieldlines become radial in the poloidal plane and that the 
flaring zone near the base of the wind is well defined. Thus 
conversely to SLIT05 where the dead zone was too ex- 
tended, we have a more realistic geometry. Figure [6] shows 
the polar radial velocity where the vertical dashed line rep- 
resents the Alfven radius, which corresponds to the border 
between application of models A and B. 

The presence of three critical points characterizes a 
different topology for this kind of wind solution (similar 
cases were already discussed in STT04). Figure [7] shows 
the profile of the temperature for this solution. It is in 



reasonably good agreement with observations, namely at 
I AU. The polytropic index is shown in Fig. [51 It shows, as 
expected, that a value around 1.2 is very well adapted to 
the solar wind, except in the low corona. However, as we 
shall discuss in next section, this first hybrid solution suf- 
fers from the same drawbacks as the previous one, despite 
its more sophisticated structure. In addition, the density 
is too low by one order of magnitude. Thus, although the 
temperature profile is low enough to be in agreement with 
observations, the mass flux at 1 AU remains too low. The 
low effective temperature does not prevent us from ob- 
taining large velocities but rather from obtaining large 
magnetic field at large distances. Thus we reconsidered 
the values of some of the parameters, in particular the 
value of the latitudinal dependence of the pressure which 
is not very well constrained from the observations, to con- 
struct another hybrid solution more fitted to the observed 
magnetic field. 

Another way of analyzing the drawback of this solution 
is to examine the convergence of the value tq /r* . It repre- 
sents by itself the convergence of the radial magnetic field 
intensity at 1 AU, Eq. (fTSll .The convergence of the values 
of the Alfven Mach number and density at this point, Eqs. 
(fT5| and (flBjl . follows as a consequence. Ulysses data at 
1 AU lead to a very high value of Rq = rQ/r^ which, in 
turn, means that the acceleration of the wind up to the 
Alfven speed should take place on a larger scale than the 
model predicts. Thus, a more satisfying solution should 
display a lower total acceleration from the surface up to 
the Alfven point. The fully radial model used (see Sect. 
[3]) is not able to produce that kind of behaviour. Some 
degree of flaring is needed in order to slow down its accel- 
eration. Hence, the only way to deal with the problem is 
to decrease the acceleration zone by decreasing the value 
of Br at 1 AU, from Eq. (jlSp. On the other hand, in model 
hybrid 1, the high density gradient near the equator (a con- 
sequence of the high values of e and 6) provides a poloidal 
pressure towards the pole, improving the polar coUima- 
tion and subsequent acceleration of the wind. Therefore, 
the equatorial gradient of the radial magnetic field must 
increase (|^| must increase), generating a higher magnetic 
pressure towards the equator which counterbalances the 
previous effect. 

This wind solution needs a value for the magnetic field 
anisotropy parameter, fj,, very different from the one ex- 
pected. It cannot easily describe the latitudinal profiles of 
the physical quantities from Ulysses data at 1 AU. Its main 
limitation (besides the temperature, to which we will re- 
fer later) is its discrepancy on the value of radial magnetic 
field. Both radial and hybrid solutions fail completely in 
reproducing the observed values of the the magnetic field 
but the hybrid solution also has a problem in reproducing 
the density at 1 AU. This lead us to construct the hybrid 
solution in a slightly different way. 
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Fig. 6. Panel (a) the polar radial velocity as a function of distance for solution hybrid 1. The vertical dashed line 
is at the transition point (Alfven point). Panel (b), a zoom view of the three critical points is displayed. The slow 
magnetosonic, the Alfven and the fast magnetosonic critical points are labelled with +, * and x, respectively. 




Fig. 7. Effective temperature as a function of distance for 
the polar axis, for solution hybrid 1. The vertical dashed 
line represents the transition point between the use of 
model A and B. 



4.3. Results for a fully converged hybrid solution - 
hybrid 2 

In Table [5] we show the input and output values of the 
most important physical quantities, for a second hybrid 
solution - hereafter referred to solution hybrid 2. Although 
we had to release some of the initial values of the param- 
eters as deduced from SLIT05, this new solution shows 
a much better agreement between the initial guesses and 
the computed values. 

For this solution, we had to change the value of k and 
of 6, although in a less dramatic way. Changing the value 
of K in the sub- Alfvenic regime is not a serious problem. 



Fig. 8. Radial profile of the polytropic index for solution 
hybrid 1. The vertical dashed line represents the transition 
point between the use of model A and B. 

First because, in the super- Alfvenic region where the value 
of K is determined for the Ulysses data, we use model B 
in which we have no control over the latitudinal depen- 
dence of the pressure. Second, fitting the value of k in this 
domain is almost irrelevant as we do not really control the 
kinetic temperature (and pressure) which is the real tem- 
perature measured by the spacecraft. This discrepancy in 
the parameter can easily be evaluated by comparing input 
and output values of the same fiow quantities. Figs.[9]and 
[Till show the geometry of the fieldlines and how it changes 
from model A to model B. A new feature of this solution 
can be seen in Fig. [TT] (a) - a zone where the radial ve- 
locity attains a local minimum, close to the Alfven point. 
Figure [12] displays the effective temperature profile along 
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the polar axis. This plot clearly shows that the kinetic 
pressure alone cannot account for the effective tempera- 
ture. Figure [T51 shows the corresponding polytropic index. 
The excess between this temperature and the observed 
one can be accounted for, as in SLIT05, assuming extra 
pressure from Alfven waves or turbulent or ram pressure 
(Fig. [T4| . We have calculated the amplitude of the ram 
velocity /magnetic field fluctuations, if the effective tem- 
perature we have calculated is assumed to be the result of 
turbulence/Alfven waves. For ram pressure we assume: 

Pram = \p^V^ , (27) 

and for the Alfven pressure we take 

-PAlfvon = -Z , (28) 

These calculations are arbitrary and a mixing of vari- 
ous processes is probably the source of the extra pressure 
that accelerate the fast wind. However it gives an order 
of magnitude of the fluctuations needed. Comparing the 
various plots a), b), c) and d) in Fig.[TH we conclude that 
Alfven waves are more appropriate to explain the acceler- 
ation in the sub-Alfvenic part where the magnetic field is 
dominant. Conversely, turbulence and fluctuations of the 
velocity may account for the acceleration in the super- 
Alfvenic region. We arrive at this conclusion only because 
the calculated fluctuations of the magnetic field in the sub- 
Alfvenic region are smaller than the calculated turbulent 
velocity field and the reverse holds in the super-Alfvenic 
part. This is the best way to minimize the amplitude of 
the fluctuations in both regions. A mixture of the two 
components is probably more realistic but this needs a 
more detailed model to interpret the role of turbulence in 
heating the flow. 

This new hybrid solution generates a field geometry 
that is continuous at the transition point (it is still not 
differentiable and kinks in the field are unavoidable) and 
shows features expected for the solar wind (Figs.OandO. 
It is also capable of reproducing almost all Ulysses data 
at 1 AU. Despite the slight difference between the values of 
the anisotropy parameters when calculated by fitting the 
data (Table [1]) and the one from the critical solution itself 
(Table [2]), most of the limitations of the previous solutions 
have been solved. Such a discrepancy should not be very 
important since the new values can easily be fitted (with 
some degree of accordance) to the observed data (see for 
instance the fit of the density in Fig. [15] and of the radial 
velocity in Fig. [16]). The ratio between the input and the 
output values for the most important physical quantities 
is very close to unity and all the continuity criteria are 
satisfied. Nevertheless, k, 5 and ^ have values departing 
from the expected ones. The value of /i is the result of the 
transition conditions stated in Eqs. (|2T]) to (|26]l and has 
been calculated accordingly. 

The values of the magnetic field intensity at 1 AU con- 
strain the value of Rq , and therefore the length of the wind 
acceleration zone (or the dead zone). Consequently, the 




Fig. 9. Fieldlines and density contours concerning the hy- 
brid solution presented in Table [2] {hybrid 2). Distances 
are given in solar radii. The three black circles represent 
the three surfaces: the slow magnetosonic separatrix and 
Alfven singularity, which are almost coincident, and a fast 
magnetosonic separatrix. 




□luatQrlBl axis 



Fig. 10. Fieldlines and density contours for solution hybrid 2, 
close to the solar surface. Distances are given in solar radii. 

physics controlling the hybrid model forced us to adapt it 
in order to obtain the required feature. Reminding that k 
and 5 characterize the anisotropy of the pressure and den- 
sity, decreasing both parameters will lead to a decrease 
of the pressure gradient towards the pole, which enables 
the wind to accelerate more slowly (from the solar sur- 
face to the Alfven point). As a consequence of accelera- 
tion continuity at the transition point, |/Lt| also diminishes 
which means that the magnetic pressure gradient towards 
the equator, outside the Alfven sphere (in the fully radial 
zone of the model), also decreases. For the dynamics of 
the radial part of the hybrid model, the wind velocity is 
expected to be higher in order to satisfy the values at 1 
AU and so it needs to accelerate the wind. This leads to 
a decrease in the magnetic pressure towards the equator 
and thus a decrease of /i. 

Of course, there is a price to pay to fit all data at 1 
AU. Thus, this hybrid model has two major drawbacks. 
The physical quantities are discontinuous at the transition 
radius except for the polar ones. This is a natural conse- 
quence of the analytical expressions that describe the flow, 
Eqs. ([T]) to and the temperature behaviour. The first 
one is solved only for values of e = 1. The second and 
more serious drawback is the very high effective tempera- 
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Fig. 11. Panel (a) the polar radial velocity as a function of distance for solution hybrid 2. The vertical dashed line is 
at the transition point (Alfven point). Panel (b), a zoom view of the three critical points: the slow magnetosonic, the 
Alfvenand the fast magnetosonic critical points, labelled with +, * and x, respectively. 




Fig. 12. Effective temperature as a function of distance 
for the polar axis, for solution hybrid 2. The vertical 
dashed line represents the transition point between the 
use of model A and B. 



turc. This can be explained only if we calculate the heat 



flux u sing a reasonable kinetic theory (jZouganelis et al 



200i) together with solving a full energy equation. This 
amounts to invoking a non thermal heating term, a dif- 
flcult task that we postpone for future work. In Fig. [T^ 
we see how the energy equation can be essential. The ab- 
sence of the an abrupt increase of the temperature very 
close to the surface in other models, such as the one pre- 
sented in SLIT05 and the one presented in Sect. [3] of the 
present work might be explained by solar surface being 
much closer to the Alfven radius and therefore the prob- 
lems had not emerged yet. Nevertheless, high tempera- 
tures are reached (for an overall behaviour of the wind 
solution) as a consequence of high values of the magnetic 
field at 1 AU and not necessarily high values of velocity 
as one may expect intuitively. 



Fig. 13. Radial profile of the polytropic index for solution 
hybrid 2. The vertical dashed line represents the transition 
point between the use of model A and B. 



5. Conclusions 

From the constrained parameters obtained after fitting 
Ulysses data (SLIT05) we were able to build different 
critical solutions for the solar wind. The first was obtained 
using a purely radial field (model B). The remaining two 
solutions where constructed as hybrid ones incorporat- 
ing an inner region where model A (with flaring stream- 
lines) was used and an outer one with model B. Thus, 
we combine the advantages of model B of reproducing 
highly adaptable functions of latitude and the advantages 
of model A of ensuring adequate flaring of the fieldlines 
to get a more realistic geometry of the overall solar wind. 
These two distinct models were coupled using well defined 
domains for each one and a suitable transition zone, the 
Alfven radius. Both model B (used by itself) and the hy- 
brid model (A and B coupled) were used to generate a 
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Fig. 15. Scaled proton density with latitude. The points 
represent hourly averaged data from the ULYSSES Swoops 
Ions experiment. The solid curve corresponds to the fit 
using model B and the bold solid curve to the solution 
generated by solution hybrid 2. 
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Fig. 16. Radial velocity with distance to the Sun. The 
points represent hourly averaged data from the ULYSSES 
Swoops Ions experiment. The curves correspond to the ra- 
dial velocity at different latitudes of the solution generated 
by solution hybrid 2. 



critical wind solution that could replicate the measured 
values of some important physical quantities, such as the 
radial velocity, the radial and toroidal magnetic fields and 
the density, at 1 AU measured by ULYSSES at solar mini- 
mum. 

In the hybrid model, we could not avoid some discon- 
tinuity out of the polar axis. This may be solved in the 
future by using the present solutions to start more realis- 
tic MHD simulations. For the first model with helicoidal 



streamlines, the computed velocity and density values at 1 
AU are in good agreement with Ulysses data. However, 
we were clearly not able to reproduce the values of the 
radial magnetic field. This is also true for the first hy- 
brid model we presented (hybrid 1), which in addition 
failed to reproduce the density at 1 AU by one order of 
magnitude. This can be understood because we tried to 
minimize the temperature along the flow. The second hy- 
brid model is able to reproduce the values of all physical 
quantities at 1 AU except for the temperature. This is too 
high even though we can invoke non thermal processes to 
explain the excess of effective pressure. It also provided 
a solution with a smooth geometry where the fieldlines 
became purely radial after the Alfvcn radius. Some con- 
cessions were made in the values of the parameters that 
rule the dynamics of both models. They are slightly differ- 
ent from the ones calculated in the method presented in 
SLIT05. Such differences can be explained by the physics 
that describes the dynamics of the flow. Even though the 
high values of the temperature at 1 AU can be explained 
(SLIT05), its behaviour close to the solar surface suggests 
the implementation of an energy equation. Since this is a 
formidable task, we postpone it for future work. However, 
we face the same difficulties as any MHD simulation. The 
main advantage of the present solutions is its simplicity. 
Moreover, they are exact solutions of the ideal MHD equa- 
tions. 
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